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Abstract 

A simple power-law growth of charged-particle multiplicities near midrap- 
idity in central Au-|-Au collisions at ^sjvjv = 56 and 130 GeV, recently mea- 
sured at RHIC, is derived. We give predictions for the central particle densi- 
ties up to y/sNN ~ 1800 GeV. A strong growth of the Au+Au densities above 
those for p{p) — p collisions is predicted. 

Charged-particle multiplicity densities near midrapidity in Au+Au collisions at 
^JSNN = 56 and 130 GeV are among the first resultsQ] from the Relativistic Heavy- 
Ion Collider (RIHC) at Brookhaven National Laboratory. The most striking feature 
of the data is the strong rise of the density (normalized to the measured number 
of participating pairs of nucleons) in going from 56 to 130 GeV, an increase of 
about 31%. As seen in Fig. 1, the rising curve through the data for Au+Au colli- 
sions is much steeper than the curve through the corresponding densities measured 
for p(p) — p collisions over the range of high c. m. collision energies, 22 to 1800 
GeV||, H, ^, 1^. The latter curve is steeper than growth as Ins. In this paper, we 
show that a simple power-law growth represents the p{p) — p central-density data 
accurately. The power is slowly increasing; the values are the same as the pow- 
ers which have been theoretically calculated on the basis of dynamical arguments 
in an early study of the average charged-particle multiplicities as a function of 
^/s. This is the curve through the p — p data points shown in Fig. 1. Using this 
power-law growth, we show that a power-law growth of the density holds for cen- 
tral Au+Au collisions, and we give explicit values for the slowly-increasing powers 
in terms of those calculated for p(j)) — p collisions Q. The power law growth fits 
the new data[D for central Au+Au collisions, and allows definite predictions for 
higher y^sjvw- These are given for ^/snn in the energy domain 200 to 1800 GeV, 
as shown by the curve through the Au+Au data in Fig. 1. We note that Pb+Pb 
dataQ at ^/snn = 17.8 GeV is also represented by the curve. It is worth recalling 
that general theoretical arguments exist for simple power-law growth of multiplici- 
ties, with a limiting power approached asymptotically, in particular the arguments 
of Polyakovjs] from the early days of approximate scale-invariance and self-similar 
processes. Also, power-law (inverse) behaviour has been derivedQ for the asymp- 
totic S-matrix amplitude whose integration over impact parameter determines the 
growth of the total p{p) — p cross section at very high energies jl^. 



The growth of the normahzed particle density near central pseudorapidity, for 
pip) ~ P collisions, is given by the form, 

^^pi^s)=AiV-s?'^^^ (1) 

where is a constant of order of unityp^, and p{s) is a small power which grows 
slowly with increasing s, approaching a value of the order of 0.1 at very high ^/s. 
This follows from an accurate representation of the average charged-particle multi- 
plicities given by the form||^, 

=A/(s)(ln.)(yi)2p(«) (2) 

Here, A is a normalization constant of order unity|^^^, f{s) is a function which 
approaches approximately unity at large s, and the (In s) factor represents the 
increasing extent of the rapidity "plateau" . Thus the essential dependence is as 
{ln.s){^/s)'^P^''\ The slowly increasing p(s) have been theoretically calculated at 12 
values of ^/s from 14 to 40,000 GeV (as given in Table I of |^). The calculation 
involves a single phenomenological parameter which controls the (slowing) growth 
of p{s) toward a limiting value of order of 0.1 (Eq. (5) in Q). The calculation is 
based upon the physical idea that collision energy which is progressively removed 
from the extreme fragmentation region sustains a power-law growth of {n{s)) with 
a slowly increasing power, i. e. an increasingly high multiplicity of particles over 
the rapidity plateau. This is consistent with the entire central p(j)) — p collision 
system becoming "blacker" over the impact-parameter plane. P|p^ We assume that 
the previously calculated powers, p{s) used in Eq. (2), are approximately the powers 
to be used in Eq. (1). Using the calculated p{s) from Table I in ||], we obtain from 
Eq. (1) the following p(yi), for A = 1.2 (with in GeV), 
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The resulting curve shown in Fig. 1 is an accurate representation of the data. Eq. (1) 
gives agreement with 7 explicit data points with a single parameter Q controlling 
the growth of the power p(s), and the normalization A. 

Given the basic power-law behavior, then for central Au-|-Au collisions, the 
growth of the normalized particle density near midrapidity can also be given ap- 
proximately by a simple power-law form, with known powers and normalization, 

1 dN 



0.5(7Vpart) dr] 

/i]^)'^(^""^ \-^{A[ ] } (4) 



{Ai^s 




^^The physical meaning of this simple normalization is about one particle per unit of rapidity 
at Vi ^ 5 GeV. [| 

^^It is noteworthy that analyses of certain characteristics of cosmic-ray air showers above 10^* 
eV suggest a more strongly growing hadronic cross section and/or "inelasticity". 
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The first quantity in brackets represents the multiphcity from each of 0.5(A^pait) 
collisions 1^ where (A'part) is the experimentally estimated average number of 
participants at ^/s'nn- The second quantity in brackets gives rise to an expanded 
multiplicity, in the approximation in which each particle from the initial multiplicity 
undergoes an additional collision^ with a typical collision energy approximated 
roughly as \/snn = {^/snn/ {n{sNN)))- This gives rise to multiplication of the 
initial multiplicity by the bracketed number which is determined by the power p, 
the value of p{s) appropriate to -n/s/vtv- Additional collisions can occur, but in the 
energy range ^/SNN from 56 to 1800 GeV discussed below, further reduction to a 
typical subsequent collision energy results in low energy. In calculating pAui\/SNN) 
from Eq. (4), we use the values of the average multiplicity {n{,/s]^)) as these are 
tabulated in Table I of We give the effective power ^ of satjv, p' = {p{snn) + 
p{snn)) and the results for pau{\/snn), 

^i]^ = 56, y ^ 0.11 ^ pAu(56) ^ 2.6 

Vi^=130, p' = 0.12 ^ pAu(130) = 3.4 

y5]^= 200, p'^S 0.125^pAu(200)?5S4 (5) 

y/JWN = 1800, p' 0.135 ^ pAu(1800) 7.1 

As a check, for Pb+Pb collisions at a relatively low ^/sWn = 17.8, p' = 0.11 — » 
/9pb(17.8) =2.0 The resulting curve in Fig. 1 represents the present the data and 
predicts the densities at higher ^/snn, from 200 to 1800 GeV. With the measured||l| 
(-/Vpart) = 330 at 56 GeV and 343 at 130 GeV, our calculated values of (dN/drj) 
are 430 and 580, respectively. The experimental numbers for {dN/drj)\,-i^n ||l| are 
respectively, 408 ± 12(stat) ± 30(syst), and 555 ± 12(stat) ± 35(syst). The predicted 
(dN/dr]) are 680 at 200 GeV (for (A^part) = 343), and 1220 (1350) at 1800 GeV 
(for (A^part) = 343(380)). In Eq. (4), the quantity {^/IWiv / {nis n n))^^ gives an 
increase in the multiplicity by about 20% at 56 GeV. This increases to ^--^30% at 
130 GeV, and to ~50% at 1800 GeV. It is useful to compare these multiplicities 
with another recent calculation jl^, which gives a low value of ~ 945 at 1800 GeV 
(i. e. pAu(1800) = 5.5). This calculation involves applying a large "correction", 
motivated by speculative dynamics, to reduce excessive multipicities from a dual 
model (already necessary at 56 GeV, i. e. note Fig. 1 in jl^). This results in a very 



strong suppression of multiplicities in the TeV range 1 12 . 

In conclusion, although scale invariance |^ is not exact, a global Q collision 
property like particle densities near midrapidity for Au+Au collisions and p{p) — 
p collisions can be directly related and quantitatively represented by a physically 



motivated, simple power-law growth with energy. This is relevant to simply 
estimating the maximum energy density 0, which is reached only for a very brief 
collision time. 



^''We approximate the multiplicity growth in the same way for hadron-hadron collisions i. e. pion- 
pion, pion-nucleon, nucleon-nucleon. At high enough energies, one might make explicit reference 
to partonic collisions. The notion of paxtonic "saturation" at very high y/s can be incorporated 
in the limiting value of the power p(s).te|| The last factor in Eq. (4) can be generalized to include 
neutrals, and a lower, effective y/s for centrally-produced particles. As an example, a resulting 
approximate numerical factor of ~ 1.5 (^Ys ^ ^) close to unity. Eq. (4) can be generalized 

to A{^/s)'^^^''^ I (1 — + xA{y/s/ (n(s)))^^ | , where the function x = for p(p) — p collisions goes 
continuously to x = 1 for central collisions of two heavy ions, at present y/s. 

^* In our calculations with all relevant v'siVJV < 50 GeV, we use an approximate p(\/sjvjv) = 
p(53) = 0.055. |q] In more detail, use of smaller powers at low energies |6| would tend to be compen- 
sated by use of higher values of effective energy for some collisions i. e. nucleon-nucleon, valence 
quarks. 

^5 Then (V«JViv/(ri(sjViv)>) is ~ 3 GeV. 

The densities as such, make no direct reference to p article momenta. 

Note that p — p data is parameterized in Q| and ll|] by very different sets of the 3 parameters 
in the form A + B{ln s) + C(ln s)^, without physical motivation. 
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Added note 



Since completion of this paper, new results have appeared from RHIC in this rapidly 
developing field. First, data on the "centrality" dependence for Au+Au collisions 
at y/SNN = 130 GeV. i. e. pAu versus (-/Vpart) appears in Fig. 4 of This can be 
approximately, but simply described by the generalized formula in footnote F3. At 
^/SNN — 130 GeV, write p = (2.25) {(1 — x) + x(1.2)(1.3)}, with the parameteriza- 
tion X = a\ii{Np/b) ^ ln(Afp/10), where Np = (iVpart). Then, 
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We assume that x remains ~ 1 for Np > 350 at y^s^vTV > 130 GeV. It is notewor- 
thy that dpAu/dNp oc 1/Np, is similar to the dpAu/dNp cx 1/Np^'^ obtained from 
the phenomenological fit pAuNp oc jYj4 i6±0.04) q^Q^-gj ^j^g above repre- 

sentation of pAuiNp) does not necessarily require a markedly increased fraction of 
"hard" processes at ^/snn = 130 GeV. In this connection, there is new data on the 
transverse energy density dEx/dri, which indicates Q that {dEr / dr]){dNch/ dr])~^ 
does not increase markedly in going down from ^/snn = 17 GeV to 130 GeV. 
The present ideas are consistent with this behavior, because the energy is going 
mainly into a marked increase in particle production (relative to p{p) — p) , not into 
a substantially larger fraction of "hard" processes. 

Finally, with respect to Fig. 1, we note that a recent measurement for Pb+Pb 
at y/sNN — 17 GeV, when scaled to Au-I-Au, gives a higher value pAu ^ 2.46. Also, 
we note that our curve for pp reaches the value of p 1.6 at y/s = 22GeV, near to 
the value measured for pp-coUisionsQ. 

We thank Dr. Mark Tannenbaum of the PHENIX collaboration at RHIC for 
stimulating correspondence. 
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Figures 




Figure 1: Measured pseudorapidity density near midrapidity normalized per mea- 
sured participant pair, for central Au+Au collisions at y^SNN = 56 and 130 GeV.||l| 
The point at 17.8 GeV is for Pb-Pb collisions. The predicted curve is calculated 
from Eq. (4). For comparision, p — p data points from ^, || are shown. This 
curve is from Eq. (1). Both curves are simple power laws with related, theoretically- 
calculated powersH, and related normalization. 



6 



